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Abstract 

This paper is concerned with the derivation of computable and guaranteed upper 
and lower bounds of the difference between the exact and the approximate solution 
of a boundary value problem for static Maxwell equations. Our analysis is based 
upon purely functional argumentation and does not attract specific properties of an 
approximation method. Therefore, the estimates derived in the paper at hand are 
applicable to any approximate solution that belongs to the corresponding energy 
space. Such estimates (also called error majorants of the functional type) have been 
derived earlier for elliptic problems \19\ [20] . 
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1 Introduction and notation 

The main goal of this paper is to derive guaranteed and computable upper and lower 
bounds of the difference between the exact solution of an electro-magneto static bound- 
ary value problem and any approximation from the corresponding energy space. We 
discuss the method with the paradigm of a prototypical electro-magneto static problem 
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in a bounded domain. The generalized formulation is given by the integral identity f l2.6p . 
We show that (as in many other problems of mathematical physics) certain transforma- 
tions of fl2.6p lead to guaranteed and fully computable majorants and minorants of the 
approximation error. However, the case considered here has special features that make 
(at some points) the derivation procedure different from, e.g., that which has been earlier 
applied to other elliptic type problems. This happens because the corresponding differ- 
ential operator has a nonzero kernel (which contains curl-free vector fields) and the set of 
trial functions in fl2.6p is restricted to a rather special subspace. For these reasons, the 
derivation of the estimates is based on Helmholtz-Weyl decompositions of vector fields, 
orthogonal projections onto subspaces, and on a certain version of a Poincare-Friedrich 
estimate for the differential operator curl. First, we show that the distance between the 
exact solution E and the approximate solution E (measured throughout the semi-norm 
generated by the operator curl) is equal to some norm of the so-called residual functional 

(cf. (13. ip ). If E satisfies the boundary condition exactly, i.e. Tt^^yE = G, then the latter 
functional vanishes if and only if curl E coincides with curl E. Lemma [10] shows that an 
error majorant can be expressed throughout a certain norm of (cf. (13. 2p ). However, in 
general, computing of this norm is hardly possible because it requires finding a supremum 
over an infinite number of vector fields. 

Theorem [T3] provides a computable form of the upper bound. The corresponding 
estimate (I3.12p shows that the error majorant is the sum of five terms, which can be 
thought of as penalties for possible violations of the relations (I2.ip -( IX4|) . It contains 
only known vector fields and global constants depending on geometrical properties of the 
domain. Moreover, it is easy to see that the upper bound vanishes if and only if E coincides 
with the exact solution E and a 'free variable' Y encompassed in the estimate coincides 
with curl E. Also, we show that the estimates derived are sharp in the sense that the 
estimates (I3.13P and (I3.14p have no irremovable gap between the left and right hand sides 
(Remark [T6l) . Finally, in Section |U we derive lower estimates of the difference between 
exact and approximate solutions. The corresponding result is presented by Theorem [5T] 
This estimate is also computable, guaranteed and sharp provided that the approximation 
exactly satisfies the prescribed boundary condition. 

Throughout this paper, we consider a bounded domain ^7 C with Lipschitz contin- 
uous boundary 7 and denote the corresponding outward unit normal vector hj n. E and 
H stand for the electric and magnetic vector fields, respectively, while e and /i denote 
positive definite, symmetric matrices with measurable, bounded coefficients that describe 
properties of the media (dielectricity and permeability, respectively). For the sake of 
brevity, matrices (matrix- valued functions) with such properties are called 'admissible'. 
We note that the corresponding inverse matrices are admissible as well. In particular, 
there exists a constant > 0, such that for a.e. x E Q 

c,\^\' < f^-\x)^ ■ ^, V^GM^ (1.1) 

By L^(fi) we denote the usual scalar L^-Hilbert space of square integrable functions over 
and by H(i7) the Hilbert space of real- valued L^- vector fields, i.e. L^(f2,]R^). For the 
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sake of simplicity we restrict our analysis to the case of real valued functions and vector 
fields. The generalization to complex valued spaces is straight forward. 

Orthogonality and the orthogonal sum with respect to the scalar product of H{Q) is 
denoted by ± and ©, respectively, i.e. if 



where A denotes Lebesgue's measure. Moreover, by ±1, (respectively (Bu) we indicate the 
orthogonality (respectively orthogonal sum) in terms of the weighted L^-scalar product 
generated by an admissible matrix u. 
Throughout the paper we will utilize the following functional spaces: 



Analogously, we define the spaces associated with the operators div and grad. Further- 
more, we introduce the spaces (containing the so-called Dirichlet and Neumann fields) 

UdA^) := H (curl°, Vl) n e'^ H (div^, Vl) 

= 1^ e I curl^ = 0, dive^ = 0, n x = o| , 

HnA^) := H (curlg, n) f] fi'^ H (div°, ^]) 

= e H(fi) I curl^ = 0, div/i^ = 0, n-/i^|^ = o| . 



Here and later on we write E G e^^ H(divQ, Q) if eE G H(divQ, Q). These are finite dimen- 
sional spaces, whose dimensions are denoted hj da and d^, respectively. In fact, these 
numbers are equal to the so-called Betti numbers of fl and depend only on topological 
properties of the domain (for a detailed presentation see [lO]). A basis of 'HD,eiS^) shall 
be given by special vector fields {Hi, . . . , -ff^^,}. 

Finally, we note that being equipped with the proper inner products all the above 
introduced functional spaces are Hilbert spaces. 

The classical formulation of the electro-magneto static problem for a given vector field 
F (driving force) and given e, reads as follows: Find a magnetic field 




H(curl,^]) 
H (curio, fi) 

H(cur^,^]) 
H(curl°,n) 




o 



C°°(i7), closure in H(curl,f2), 
H(curr,fi) n H (curio, fi). 



H G H(curl,fi) n jjL 



lH(div°,fi)^•H;v,M(^]) 



and a corresponding electric field 



EeH (curl°, divo e, ±e, f^) := H (curl°, n e'^ H (divg, fi) n HdA^) 
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such that in Q 

curl H = F, curl E = fiH. 
In other words, the problem is to find vector fields H G H(curl, Q) fl fi^^ H(div, Q) and 

E G H(curl,dive,fi) := H(curl, fi) n H(div, fi), 

such that 

curl H = F, curl E = fiH in Q, 

div fiH = 0, div eE = in fi, 

n ■ fiH\^ = 0, n X = on 7, 

where the homogeneous boundary conditions are to be understood in the weak sense. 

This coupled problem is equivalent to an electro-magneto static Maxwell problem in 
second order form, which in classical terms reads as follows: Find an electric field E in 
H(curr, divoe, J-e, such that fi~^cm\E belongs to H(curl,f2) and 

curl curl E = F 

holds in Q, i.e. find E G H (curl, dive, fi), such that yU'^curli? G H(curl,n) and 

curl /i"^ curl E = F in Vt, (1.2) 

dweE = inn, (1.3) 

n X E\^ = on 7, (1.4) 

eE ± HdA^)- (1-5) 

Once E has been found, the magnetic field is given by H := curlE'. 
We note that the problem 

curl curl E + k^E = F in fl, 

n X E\^ = on 7 

with positive k, was considered in [2J in the context of functional type a posteriori error 
estimates. From the mathematical point of view, this problem is much simpler as the 
problem fll.2p - fll.5l) since the zero order term makes the underlying operator positive 
definite. 



2 Variational formulation and solution theory 

Henceforth, we consider fll.2p - fll.5p assuming that the boundary condition on 7 may be in- 
homogeneous (physically, such a condition is motivated by the presence of electric currents 
on the boundary). Hence, we intend to discuss the following prototypical electro- magneto 
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static Maxwell problem in second order form: Find an electric field E in H(curl, dive, Q), 
such that 

curl curl E = F 

dweE = 
n X E\^ = G 

eE ±nD,em, 

i.e., find E in 

H (curl, divo e, X^, ^) ■= H (curl, fi) n e'^ H (divp, fi) n UdA^)^" 

satisfying (12. ip and (12. 3p . There are at least two methods to prove existence of the 
solution. One is based upon Helmholtz-Weyl decompositions (see, e.g. P Ull [151 HZl 
[Tot [TT|). The second method consists of introducing and studying a suitable generalized 
statement of the problem (I2.ip -f l?!^ . In this paper, we use the second method because it 
provides a natural way of deriving error estimates. Both methods are based on Poincare- 
Friedrich estimates, see Remark IH and (if it is needed) exploit suitable extension operators 
for the boundary data. On this way, we also need a certain version of the Poincare- 
Friedrich estimate, namely 

i^if^ < Cpicurl^l^ G H(curr,divoe,±e,f^)- (2.5) 

Of course, there exist more general variants of Poincare-Friedrich's estimate (12. 5 p for 
vector fields. Here, we refer to Remark [HI 

Now, let E^ be some vector field in H (curl, divg e, -L^, fl) satisfying the boundary con- 
dition (12. 3 p in the generalized sense, i.e., E — E^ & H(curr,r2). The generalized solution 

E (curr, divo e, ±£, fi) + C H (curl, divo e, X^, fi) 

of (12. ip - (12.41) is then defined by the relation 

{l2~\m\ E, curl W) ^ = {F,W)^ VW^ G H(curr, divo e, ±„ fi). (2.6) 

If F G H(f2) then by the Cauchy-Scharz inequality the right hand side of (12. 6p is a 
linear and continuous functional over H(curr, divo^, J-e, fi). By (12.50 the left hand side 
of (12. 6p is a strongly coercive bilinear form over H (curl", divg e, ±£, Q). Thus, under these 
assumptions the problem (12.60 is uniquely solvable in H(curr, divo^, J-e, fl) + E^ by Lax- 
Milgram's theorem. 

First, we note some Helmholtz-Weyl decompositions of H(f2), i.e. decompositions into 
solenoidal and curl-free fields, which will be used frequently throughout our analysis. 

Lemma 1 H(fi) can be decomposed as 

H{n)=eH (curio, n) curl H (curl, fi) 
= £:gradH(grad°,fi) ©^-i H(divo,fi) 

= £gradH(grad°,n) ©^-i sHD^eiP) ®e~^ curl H (curl, fi) (2.7) 



in (2.1) 
in Vt, (2.2) 
on 7, (2.3) 
(2.4) 
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and 



H (^]) = H (curl°, n) ^-^curl H (curl, Q) 
= gradH(grad°,l]) e'^ H (divQ, n) 

= gradH(grad°,fi) noA^) ^"^curl H(curl, fi), (2.8) 

o 

where all closures are taken in and H(grad°,r2) = H^(r2). Moreover, 

curl H (curl, fi) = H (divo, ±, fi) := H (divp, n) n HdA^)^- 

Remark 2 Let us denote the e-orthogonal projection onto e^^curl H (curl, fi) in fl2.8p by 
71. Then we have for all $ G H(curl, Q) 

n_^7r$ = rt^^$, curl 7r$ = curl $ (2.9) 

and for all eH{Q) 

div£7r^ = 0, 67i'^±nD,e{^), curl(l - tt)^ = 0, n,^(l - tt)^ = 0. 

The latter line can be written in a more compact and precise way as 

7rH{Q) =e-^curlH(curl,fi) = H (divo £, ±e, f^), {I - = H{cmll,Q). 

Remark 3 Note that by (12.11) F must be solenoidal and perpendicular in H{Q) to 
H(curlQ,fi). Using the Helmholtz-Weyl decomposition (12. 7p we decompose the vector field 
F e V\{VL), i.e. F = eFd + e-Fgrad + -Fcuri- Then, for any W (curf, divo e, ±£, Vt) we 
compute {F, W)q = (-Fcuri, W)n- Hence, the functional on the right hand side of (12.61) can 
not distinguish between F and the projection Fcuri- 

The following theorem states the main existence result. 

Theorem 4 Let F G H(divo, -L, fl) and let G H(curl, divoe, J-e, Vt) satisfy the bound- 
ary condition (12. 3p . Then the boundary value problem (12. ip - (12.40 is uniquely weakly solv- 
able in H (curr, divo e, -L^, Q) + E^. The solution operator is continuous. 

Remark 5 The kernel of (12. II) - (12. 3p equals T-Lo.ei^)- ^^^?/ have to show curli? = 
but this follows immediately since E G H(curr,fi) and thus, 

= <^curl fi~^ curl E, i?)^ = (yU~^ curl E, curl i?)^ . 

Remark 6 The boundary data G and its extension E^ can be described in more detail. 
Since the papers IH |^ and the more general paper of Week we know that even 
for Lipschitz domains, where the non scalar trace business is a challenging task, there 
exist a bounded linear tangential trace operator Tt^^ and a corresponding bounded linear 
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tangential extension operator ft^^ (right inverse) mapping H(curl,f2) to special tangential 
vector fields on the boundary, i.e. 

H7^/'(curl„7) := |^ G \ curl.V^ G H-V2(^)| , 

and vice verse. Here, curls denotes the surface curl. Using the Helmholtz-Weyl decompo- 
sition fl2.8p we even get an improved extension operator. We have 

— 1/2 

Tt^^ : H(curl,fi) ^ (curls, 7), 

— 1/2 

n,^ : Hi (curls,7) H(curl,divoe, ±£,11). 

Applied to smooth vector fields we have Tt^^ = nx -1^. Now, we may specify the boundary 
data G E (curls, 7) and the extension := ft^^G G H(curl, divo e, -L^, fi) as well as 
our variational formulation for E = E° + ff^^G: Find E° E H (curl", divg e, J.^, Q), such 
that 

b{E°, W) := {fi-^ curl E°, curl W)^ = (F, W)^ - (//"^ curl fi,^G, curl W)^ =: i{W) 
holds for all W G H(curr, divo e, -Le,^)- 

Remark 7 Henceforth, we assume that G is given by a tangential trace of some vector 
field G G H(curl,n). 

Remark 8 More general variants of the Poincare-Friedrich estimate for vector fields 
(12. 5p are known. For instance, we have 

mn < ||curlvl>||^ + Wdwen^ + K^n^-m^^^^^^^^^ + 5^ |(£Vl/,i7„)^ 

\ n=l 

which holds for all \1' G H (curl, dive:, 1]). This estimate may be proved by an indirect 
argument using a Maxwell compact embedding property' ofQ, which holds true not only 
for Lipschitz domains, but also, if the homogeneous boundary condition is considered, 
for more irregular domains (cone properties), see fT^ . For inhomogeneous boundary 
conditions the Lipschitz assumption can not be weakened. Actually, it is just the continuity 
of the solution operator of the corresponding electro static boundary value problem, see 

3 Upper bounds for the deviation from the exact so- 
lution 

Let E be an approximation of G H (curl", divo e, _Le, Q) + E^ C H (curl, divo e, _Le, Q). We 
assume that E belongs to H (curl, dive:, fi), which means that, in general, the boundary 
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condition, the divergence-free condition, and the orthogonahty to the Dirichlet fields might 
be violated, i.e. the approximation field may be such that 

Tt^^E ^ G, diYsE ^ 0, (^eE, ^ for some H G UdA^)- 

Moreover, for the subsequent analysis and then also for the numerical application, which 
is even more important, it is sufficient to assume just E G H(curl, Q). 

Our goal is to obtain upper bounds for the difference between curli? and curli? in 
terms of the weighted norm 

1*1-. lk-"^'p|L = (^'-''p.'p>f ■ 

First, we use (12. 6p and get for all W G H(curr, divoe, J-e, VL) 

^i-^c\xi\{E - E),c\xi\w\ = {F,W)^- U-^c\xi\E,cm\w\ =: ^^(ly), (3.1) 

where ^^ is a linear and continuous functional over H (curl", divo e, J-e, i^) as well as over 
H(curl,fi). Furthermore, does not depend on the exact solution E. 

Remark 9 Obviously, vanishes ifcmlE = curli?. Furthermore, if E satisfies the 
boundary condition exactly, i.e. Tt^^E = G, then = Q if and only if cmlE = curli^ (or 
what is equivalent if and only if E = ttE). This holds by the following argument using the 
Helmholtz- Weyl decomposition: If Tt^^E = G then E — ttE G H (curl", divo e, ±e, fi) . Thus 
curl(ii^ — 7tE) = by = 0. But then E — ttE is a Dirichlet field and hence must vanish 
by orthogonality. Finally cutIttE = curli?. 

The second step is based upon the following result: 

Lemma 10 Let E G H(curl, divo J-e, f2) be the exact solution and E G H (curl, 17) be 
an approximation. Furthermore, let be as above and let q > exist, such that 

jj-^ cm\{E - E), curl W^^ = l^iW) < a fcurl Wj^^^^^ 

holds for all W E H (curl", divg e, ±e, 17) . Then 

cml(E-E) < Q + 2||curlT|l i„ (3.2) 

fi-\n ^ 

holds for all T G H (curl, 17), for which the tangential trace coincides with the tangential 
trace of E — E, i.e. G — Tt^^E, on the boundary 7. // additionally Tt^^E = G then 

curl(E -E) < Q. (3.3) 
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Proof We use the Helmholtz-Weyl decomposition fl2.8p and the projection vr from Remark 
m We consider a vector field T G H(curl, Q) with Tt^^T = G — Tt^^E and define the vector 
field 

W := E - tt{T + E) = E - E + {1 - tt)E - 7iT e H (curf, divo e, ±e, fi), 

which holds by (12. 9p . Hence, curl = curl(i? — E) — curlT. Using Cauchy-Schwarz' 
inequahty we obtain 



II curl ly I 



cni\{E - E), curl W^^ - {/j,-^ curlT, curl W)^ 
< + ||curlT||^_i_j^) ||curliy||^_i^^ 



and thus fcurl TV ||^_i ^ < Q + ||curlT||^_i ^. By the triangle inequality we get (13. 2p . (13.31) 
is trivial setting T := 0. □ 

Using the trace and extension operators from Remark E] we obtain the following result: 
Corollary 11 Let the assumptions of Lemma UU be satisfied. Then 



curl(E - E) 



< Q + 2 



< Q + 2c^ 



G - Tt^^E 



(3.4) 



Ht-l/^(curl,,7) 



Here > zs the constant in the inequality 



||curlfi,^?/'||^_i^^ < 



H-/^(curi.,7) V^/^e H, '/'(curl3,7). 



(3.5) 



Proof Setting T := ft^^{G — Tt^^E) in (13. 2p and using (13.51) proves (13.41) . We note that 
(13. 3p follows directly from the corollary as well. □ 

Lemma [TUl and Corollary [TT] imply the following result. 

Theorem 12 Let E,E be as in LemmalTR Then 

cm\{E - E) 



< 



\\F - curlFfo + 



curlE 



+ 2c., 



G - Tt.,E 



(3.6) 



H, (curl,, 7) 



where Y is an arbitrary vector field in H(curl,f2). 
Proof For any Y G H(curl, Q) and any W G H(curr, Q) we have 

- (curl Y, W)^ + {Y, curl W)^ = 0. 



(3.7) 
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Combining f l3.ip and f l3.7p . we obtain for all G H (curr, divg e, ±£, Q) 
(^;U~^curl(E - E),cutIW^ 
= (F - curl y, W)^ + (y - curl E 



I n 

,-1 curl ) 



/ n 



(3.8) 



By Cauchy-Schwarz' inequality, Poincare-Friedrich's estimate fl2.5p and (11. ip we estimate 
the right hand side isiW) of 



\{F - cmlY,W)^\ < ||F-curir||^||l^||f^ < Cp ||F - curl || curl 



<^||F_curlF||J|curlW^||^_,_^, 



Y — fi ^ curl E, curl W 



< 



Y - fi-^cmlE 



II curl I 



Now, Lemma [TD] and Corollary [TTbomplete the proof. 



(3.9) 

(3.10) 
□ 



We remark that the latter estimate is unable to measure adequately the deviation 
of the divergence of eE to (this is obvious since eE even does not need to have any 
divergence). On the other hand, even if div eE ^ then the semi-norm ||curl ■ ||^_i ^ could 

not feel the lack of the constraint diveii^ = 0. The same holds true for the deviation of 
eE from the orthogonality to the Dirichlet fields. However, it is not difficult to transform 
the estimate into a form, in which the estimate is represented in terms of the semi-norm 



do 



II ^11 



curl^||^-i,f^ + ||div£^||f^ + ^ |(e^, Hn)i 



(3.11) 



n=l 



on H(curl, dive, Q), which obviously is a norm on 

H(curr,dive,fi) := H(curr, fi) n e"^ H(div, fi). 

Remark 13 These two facts can be seen applying the Helmholtz-Weyl decomposition 
(12. 8p and the projection 7T . In particular, by (12.91) replacing E bynE in TheoremlTB would 
change nothing. In other words, the part of E containing the eventually non vanishing 
divergence term (1 — 7r)E can be added to any term in (13. 6p without changing anything. 
To get the 'full' norm (13. lip we just add the terms 



and the sum of 



dive{E - E) 



eiE-E),H„ 



diveE 



(eE,H, 
\ / n 



dived - T^)E 



e{l - 7r)E, H„ 



n 



Of course, the terms in the first equalities make sense for eE G H(div, Q) only. 
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Theorem 14 Let E, E be as in LemmaUU and additionally E G H(curl, dive, Q). Then 



E-E 



<M+{E-Y) :-- 



IF - curl F II 



Y - fi'^cmlEl 



G - r.^E 



,-1/2 



(curls, 7) 



+ 



diveE 



do 



n=l 



(3.12) 



holds for any Y G H(curl,f2). If E — E even belongs to H(curr, dive, i.e. if the 
approximation E satisfies the boundary condition exactly, then ||| ■ is a norm for E — E 
and we have for a// F G H (curl, Q) 



E-E 



<M^{E-Y) 



\\F - curl F II 



Y - /i"^ curl E 



+ 



div sE 



do 



n=l 



(3.13) 



Remark 15 If E satisfies the prescribed boundary condition and eE is solenoidal and 
perpendicular to Dirichlet fields, then (13. 6p or (I3.12p . (13.131) imply for all Y G H(curl,f2) 



E-E 



curl(E - E) 



<M+{E-Y) 



||F - curlFllo + 



F-/i"^curlF 



(3.14) 



and the left hand side is a norm for E — E . The estimates (I3.6p -( 1XT^ show that deviations 
from exact solutions contain weighted residuals of basic relations with weights given by 
constants in the corresponding embedding inequalities. These are typical features of the 
so-called functional a posteriori error estimates. 

Remark 16 We see that M+{E; Y) = 0, if and only if E := E and Y := cnilE (by 
Lemma fTPI we then have Y G H(curl,fi)^. Moreover, we note that (I3.13P and (I3.14p are 
sharp, which easily can be seen by setting Y := fi'^cuilE G H(curl,fi). In other words, 
if E E H (curl, dive, f2) satisfies the boundary condition exactly then 



E-E 



inf M:(E:Y). 

yGH(curl,n) 



Remark 17 In Theorem\T^ and TheoremYT^ we can replace the boundary term on the 



right hand side by 2||curlT|| 



or 2 



using Lemma [73 and 



curl ft^y{G — Tt^-yE) 

Corollary [771 Especially for numerical applications the first choice is recommendable. 
Hence, we may assume that G is always given by a tangential trace of some vector field 
G G H(curl,fi), i.e. Tt^^G = G. Then T := G — E and we do not have to know the 



constant c^. 



Remark 18 // the domain is 'simple' in terms of a vanishing second Betti number, i.e. 
there are no 'handles', then there exist no Dirichtlet fields. Thus, for instance, in Theorem 



14 the last summand in the respective estimates does not occur. 
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4 Lower bounds for the error 



Now, we proceed to derive computable lower bounds of the error. First, we present 
the following subsidiary result: 

Lemma 19 If E satisfies fl2.6p then fx^^cmlE G H(curl,f2) and curl curl ii^ = F. 
Proof We need to show that 

(/i-^ curl curl = V$ G C°°(r]). (4.1) 

Using TT from Remark [21 we obtain W = 7r$ G H (curl", divo £, J-e, fi) provided that 
$ G C'^{n). Thus, by O and the fact that curl(l - 7r)$ = 0, we get 

curl curl <l>)^ = (/i-i curl curl 7r<l>)^ = (F,7r<l>)^ V<l> G C°°(^]). (4.2) 



Since F G H(divo, -L, fi) = curl H(curl, f2), we get (by approximation) {F, vr$)j^ = {F, 
and (14.11) follows. To be more precise, we select Fn G H(curl,fi), for which (curlF„)„gN 
converges in H(i7) to F, using 7r$ G H(curr,f2) and curl(l — 7r)$ = 0. Then 



(curlF„,7r<l>)^ = curl7r<l>)^ = (F„,curl$)^ = (curlF„,<l>)^ V<l> G C'^(fi). 



□ 



Lemma Uni implies 

Remark 20 Let E G H(curl,r2) and some F be given. Then the following three asser- 
tions are equivalent: 

(i) /i^^ curl E E \-\ (curl, fl) and curl fi~^ curl E = F. 
(ii) F eH{n) and 

(^"^ curl curl = V$ G H(cur^,^]). 

(iii) Fg H(divo,±,fi) and 

curl curl = {F,^)^ V<l> G H(curr, divoe, fi)- 

Theorem 21 Lei G H(curl,fi) be an approximation. Then 



> snpA'F{E;W), 



cm\{E - E) 

M.{E; W) ■.= 2{F,W)^- (/i"^ curl(2E + W^), curl 



where 



and the supremum is taken over H(curr,fi). This estimate is sharp if E — E belongs to 
the latter space, i. e. if the approximation E satisfies the boundary condition exactly. 
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Proof We start with the obvious identity 



curl(E - E) 



yGH(n) 



sup (2(n-^cm\{E-E),Y) - \\Y 



Thus, for all W G H(curl, fl) we obtain the estimate 



curl(E - E) 



> 2 curl(E - E), curl w) - fcurl Wf.^ 



2(/i"^curlE,curliy)^ - cm\{2E + W),cm\W 



Clearly, this estimate is sharp because we can always put W = E—E. However, to exclude 
the unknown exact solution E from the right hand side we need G H (curl", divo e, _Le, Q). 
Then, by ([22]) 

(/i-^ curl E, curl W)^ = (F, W)^ (4.3) 
and by Lemma OS] (gSD even holds for all W e H{cm\°,n). Thus, for all W G H(curr,fi) 



curl{E-E) >M4E;W). 

Obviously, this lower bound is sharp if we can set W = E — E E H(curr, fl). 

The following result is trivial: 
Corollary 22 Let E G H (curl, dive, fi) be an approximation. Then 

do 



(4.4) 

□ 



E-E 



> sup M.{E;W) + 

^ VKeH(curl°,n) 



diveE 



J2{eE,Hn 



n=l 



Remark 23 In general, the lower bound is not sharp because if E — E ^ H(curr,fi) 
then we can not put W = E — E anymore. In fact, with fi~^ curl E G H(curl,f2) and 
curl curlii^ = F by LemmalWi we get for all W G H(curl, Q) 



[jj-^ curl E, curl W)^ = (F, W)^ + {n,^!^-^ curl E, n^/W)^ 



(4.5) 



Here, we introduced a second tangential trace ft,^, called the normal trace in terms of differ- 
ential forms, mapping again H(curl,fi) to special tangential vector fields on the boundary 
as well, i.e. 

H7^/^(div„7) := G H;^'\^) I divs^ G H-V2(7)} , 

where divs denotes the surface divergence. In and, more general, in I2&^ . see also 

[T^ Lemma 3.7, q = 1], it has been pointed out that even for Lipschitz domains the 
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integration by parts formula fl4.5p remains valid in some sophisticated sense. For smooth 



vector fields we have Tt, 



>7 



-n X (n X 



l7- 



Hence, we obtain the estimate 



curl(E -E) > M_ {E; W) + 2 curl E, n^^W)^ 



for all W G H(curl,fi), which is sharp and coincides with (14. 4p if W ^ H(curr,f2). 
But the unknown exact solution E still appears on the right-hand side, i.e. the second 
tangential trace of fi^^cnilE on 7. Furthermore, if the term {ft^.y n"^ curl E,t^W)^ is 
positive then (14.31) can not be sharp. 
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